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REPORT  SUMMARY 


This  is  the  Final  Technical  Report  describing  the  work  done  and 
the  accomplishments  of  the  research  effort  for  ONR-ARPA  titled  "Composite 
Design  Synthesis". 

The  research  carried  out  deals  with  the  problem  of  design  syn¬ 
thesis  in  heterogeneous  elasticity.  Design  synthesis  is  defined  as  the 
achievement  of  a  desired  design  criterion,  i.e.,  stress  distribution, 
strength-to-weight  ratio,  etc.,  by  preselecting  a  stress  or  displacement 
pattern  in  a  stretched  plate  and  then  determining  the  variation  of  the 
elastic  moduli  that  is  required  to  permit  the  desired  effects.  This  accom¬ 
plishment  requires  the  solution  of  the  governing  equations  of  elasticity, 
particularly  the  compatibility  equation,  in  terms  of  preselected  stress 
fields  in  the  body  of  the.  plate  for  unknown  material  properties  which  are 
spatial  functions. 

During  this  work,  solutions  to  the  moduli  variation  problem  for 
annular  disks  have  been  achieved  for  two  stress  criteria;  constant  hoop 
stress;  and  constant  ln-plane  shear  stress.  The  disk  may  be  rotating  and 
have  boundary  traction.  A  computer  program,  D0M0V1,  was  developed  to  carry 
out  these  solutions.  Attempts  to  solve  the  moduli  variation  problem  for 
a  hole  in  an  infinite  plate  (two  dimensional),  subject  to  certain  stress 
distributions  w»re  unsuccessful.  However,  great  insight  was  gained  into 
this  problem  for  future  work. 

Basically,  this  work  has  shown  that  the  concept  of  design  syn¬ 
thesis,  as  defined  here,  is  a  workable  discipline  and  in  the  case  of 
rotating  aunular  disks  and  pressurized  thick-vall  cylinders,  can  be 
applied  utilizing  the  present  state-of-the-art  fabrication  technology, 
but  that  its  application  to  complex  problems  requires  additional  work. 
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INTRODUCTION 

It  has  long  been  recognized  that  structural  elements  composed  of 
composite  materials,  such  as  glass,  boron,  carbon,  or  other  filaments,  embedded 
in  a  suitable  matrix,  such  as  epoxy  or  polyester,  offer  outstanding  strength- 
to-weight  ratios.  The  potential  of  such  materials  is  considered  so  great  that 
material  scientists  and  engineers  believe  that  they  will  form  the  bulk  of  the 
structural  materials  of  the  future.  Though  the  application  of  such  materials 
has  been  a  growing  part  of  the  state  of  the  art  for  structural  components, 
particularly  in  the  aerospace  industry,  the  translation  of  the  concept  of 
fibrous  composites  into  a  primary  load  carrying  structure  has  been  and  remains 
a  challenging  process. 

The  present  and  growing  use  c f  structural  elements  fabricated  from 
composite  materials  creates  the  need  for  the  development  of  a  rational  analytic 
design  basis,  which  to  a  great  extent  is  presently  non-existent.  This  is  not 
to  be  construed  as  meaning  that  little  or  no  research  on  composite  materials 
has  been  carried  out.  On  the  contrary,  a  large  amount  of  literature  has  been 
generated  dealing  with  both  the  determination  of  the  mechanical  properties  of 
these  materials  and  the  analysis  of  specific  structures  fabricated  from  them. 

In  general,  from  the  microscopic  viewpoint,  research  on  the  mechani¬ 
cal  properties  of  composites  has  dealt  with  the  determination  of  such  properties 
for  materials  having  given  component  elements  ordered  in  fixed  spatial 
relationships.  The  spatial  relation  in  these  cases  might  have  a  high  degree 
of  symmetry,  as  in  long  or  continuous  filament  composites,  or  a  completely 
random  or  homogeneously  disordered  array  as  usually  employed  in  short  carbon 
or  boron  fiber  composites.  Such  research  has  been  directed  towards  the 
creation  of  analytic  or  experimental  methods  of  determining  the  mechanical 
properties  of  composite  structures  in  terms  of  the  known  properties  of  the 
composites'  components,  characteristic  of  fMo  approach  are  the  works  of 
Sayers  and  Hanley  [l]*,  Chen  and  Cheng  [2],  Hill  [3],  and  Gaonkar  [A],  among 
others. 

In  the  analysis  of  structures  composed  of  such  composites,  the 
material  has  generally  been  treated  as  exhibiting  gross,  homogeneous, 

«j|f 

Numbers  in  brackets  are  references  found  at  the  end  of  this  report. 
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isotropic  or  anisotropic  mechanical  properties.  These  gross  properties,  when 
entered  into  the  constitutive  equations  defining  the  material,  have  allowed 
analyses  of  such  st  uctures  through  the  classical  methods  of  the  theory  of 
elasticity,  plates  and  shells,  vibration,  and  others.  Along  these  lines,  the 
concept  of  the  "unidirectior.il  lamina"  was  introduced  and  utilized  as  the 
"fundamental  unit  of  material"  in  design  and  analysis.  This  procedure 
employs  test  data  obtained  from  a  unidirectional  lamina  as  the  basis  for  the 
design  of  laminated  components  and  structures.  Exemplary  work  done  along 
these  lines  has  been  carried  out  by  Dong  [5],  Tsai  [6,7],  Tsai  and  Azzi  [8], 
Whitney  [9],  and  Whitney  and  Leissa  [10,11],  also  among  others. 

All  of  these  analyses  have  dealt  with  materials  and  structures  that 
have  predetermined  mechanical  and  behavioral  characteristics.  That  is,  once 
the  geometric  array  and  the  constituents  of  the  composite  are  prescribed, 
then  the  mechanical  properties  of  the  material  and  the  response  characteristics 
of  a  structure  fabricated  from  such  a  material  have  been  inherently  established. 
Ani lysis  will  merely  determine  what  these  properties  and  response  charac¬ 
teristics  are. 

When  dealing  with  composite  materials,  the  analytical  procedures 
discussed  above  appear  to  be  highly  inefficient  in  many  applications.  The 
designe*.  of  fibrous  composite  structures  is  presented  with  numerous  degrees 
of  freedom  and  an  opportunity  to  exercise  ingenuity  totally  unavailable  to  him 
with  conventional  materials.  Composites,  whether  filamentary,  fibrous,  or 
sintered  or  fused  metallics,  are  capable  of  being  tailored  to  meet  specific 
requirements.  When  considering  specific  structural  applications  for  such 
materials,  it  would  be  logical  to  assume  that  a  structure  could  be  optimized, 
depending,  of  course,  upon  the  optimization  criterion,  by  varying  the  mechani¬ 
cal  properties  of  the  material  throughout  the  structure.  Further  it  would  be 
logical  to  bypass  analysis  completely  and  define  this  now  nonhomogeneous 
structure  by  some  means  of  design  synthesis.  Admittedly,  the  creation  of  a 
design  synthesis  procedure  to  adequately  handle  most  problems  in  structural 
design  is  quite  difficult.  However,  the  concept  of  design  synthesis  tc  deter¬ 
mine  the  variation  of  the  mechanical  properties  of  a  material  within  a  structure 
so  as  to  achieve  a  desired  stress  or  deformation  pattern  ir.  that  structure  is 
one  capable  of  being  developed. 
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GENERAL  DISCUSSION 


Consider  the  following  question: 

Given  a  plane  elastic  body  with  known  boundary  tractions  and/or 
displacements,  can  the  mechanical  properties  of  the  contiruum  be  described 
such  that  an  "arbitrary"  stress  distribution  within  the  body  is  met? 

The  term  "arbitrary"  is  to  be  understood  as  defining  a  family  of 
stress  distributions  that  are  preselected  but  still  conform  to  equilibrium 
requirements  and  boundary  conditions.  To  answer  this  question  we  start  by 
making  the  following  two  basic  assumptions: 

(1)  the  classical  equations  of  linear  elasticity  are  valid  in  this 
application,  and 

(2)  the  mechanical  properties  of  the  continuum  can  be  expressed  as 
spatial  functions. 

Following  from  these  assumptions  the  well  known  governing  relations  of  gener¬ 
alized  plane  stress,  given  in  rectangular  coordinates  are  as  follows: 


Equilibrium  Equations 


+  X  «=  0 


+  Y  «  0 


> 


Strain-Displacement  Equations 

du  dv 

ex  "  Sx  ’  ey  "  By  * 


av.au 
ax  ay  » 


Compatibility  Equation 
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xv 
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Assuming  tnat  the  continuum  exhibits  orthotronic  material  proper¬ 
ties  and  neglecting  time  and  strain  rate  effects,  the  constitutive  equations 
can  be  expressed  as  generalized  Hooke's  Law  as 


where  and  are  the  coefficients  of  thermal  expansion  and  T  represents 
the  temperature  difference  distribution.  The  relations  defined  by  Equations 
(4)  are  valid  when  the  axes  of  the  material  properties  of  the  continuum  are 
coincident  with  the  axes  selected  for  the  differential  equations  of  the 
problem.  If  the  axes  are  not  coincident,  except  for  the  z-axes,  and  the 
other  two  axes  of  the  material  properties  are  rotated  about  the  z-axis  through 
some  angle,  «(,  in  relation  to  the  geometric  axes,  then  more  complicated 
relations  between  the  stresses,  temperature  and  strains  are  developed. 
Lekhnitskii  [12]  presents  these  relationships  in  some  detail.  For  th'  gener¬ 
alized  plane  stress  case  in  point  the  material  coefficients  are  related  to  the 
two  axis  system  by 


where  m  *  cos  and  n  ■  sin  i. 

No  generality  will  be  lost  by  continuing  with  the  constitutive 
equation  as  given  by  Equations  (4).  Substituting  Equations  (4)  into  (3), 
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Making  the  appropriate  substitutions  into  Equation  (6)  yields 

{■ 
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and,  of  course,  the  equilibrium  equations  are  met  exactly. 

Notice  that  in  Equations  (4),  (5),  and  (9),  a^2  has  not  been 
equated  to  a2^.  Norm' lly,  under  the  limit  of  small  displacement  theory 
dealing  with  linearly  elastic  materials  which  are  conservative,  the  material 
coefficient  matrix,  defined  in  Equation  (4),  wou"d  be  symmetric  and  a^2  would 
equal  a21<  However,  some  recent  work  by  Bert  and  Guess  tl3 ] ,  among  others, 
shows  that  there  exists  experimentally  derived  data  which  indicate  that  for 
some  types  of  composite  materials  exhibiting  orthotropic  properties  the 
material  coefficient  matrix  is  not  symmetric  and  a^2  is  not  equal  to  a^. 

To  limit  the  growing  complexity  of  this  work,  the  material  coefficient  matrix 
will  be  taken  as  symmetric,  at  least  for  the  initial  phase  of  this  effort. 

Equation  (9)  can  also  be  expressed  in  polar  coordinates  as 


(Note:  Here,  the  a..,  are  in  polar  coordinates  and  are 

lj  s 

not  equal  to  the  a..,  for  rectangular  coordinates.) 

lj  8 


follows: 


9 


with  Y  being  the  stress  function  defined  by 


a  .  1 21  +  i.  ix  +  v 

r  r28e2  ■ 


a  "=  9  +  VQ 

9  dr  9 


ar0 


(ID 


Equations  (11)  satisfy  the  equilibrium  equations  formulated  in  polar  coor¬ 
dinates  for  plane  stress  which  are 


da  ,  da,  a  -  a, 
_ I  +  I  r9.  ,  _r _ g. 

dr  r  90  r 


+  R  =  0 


acrr6  ,  1  *°9 
dr  r  d0 


+  2 


0  . 


(12) 


In  order  for  Equations  (11)  to  meet  Equation  (12)  exactly  the  body  functions 

V  and  V  must  be  defined  as 
T  0 


v 


0 


dr 


(rVr)  ■=  R 


90 


0 


(13) 


thus  putting  a  rather  narrow  interpretation  on  the  body  forces. 

In  the  classical  approach  to  the  problem  of  orthotropic  plane 
elasticity,  the  coefficients  a^  are  either  constants,  as  in  the  case  of  the 
homogeneous  condition,  or  as  in  some  rare  cases,  are  special  functions  of 
position.  In  the  first  case  all  the  terms  within  the  second  set  of  large 
braces,  |  j-,  in  Equations  (9)  and  (10)  are  zero  leaving  the  remainder  of 
these  two  equations  in  the  form  of  the  well-known,  homogeneous,  orthotropic. 


v 


* 


10 


compatibility  equations  in  terms  of  the  stress  function  Y.  In  the  secord 
case  all  the  partial  derivatives  of  the  a^'s  which  appear  in  these  second 
sets  of  braces  are  capable  of  being  evaluated,  resulting  in  extremely  compli¬ 
cated  fourth  order  partial  differential  equations  with  variable  coefficients. 
Proceeding  along  classical  lines,  these  equations  must  be  solved  for  Y  which 
contains  arbitrary  constants  of  integration.  These  constants  are  then  deter¬ 
mined  by  evaluating  Y  in  terras  of  the  stresses  on  the  boundaries. 

Suppose  it  is  assumed  that  the  material  coefficients,  the  a's. 

ij  » 

are  unknown  but  that  the  stress  function  Y  is  a  fully  defined  function  of  the 
spatial  coordinates.  That  is,  the  stresses  throughout  the  body  as  well  as  on 
the  boundary  are  known.  In  such  a  case,  Equations  (9)  and  (10)  reduce  to 
second  order  partial  differential  equations  with  variable  coefficients,  in 
terms  of  the  a^'s.  The  solution  of  these  equations  and  the  resultant  deter¬ 
mination  of  the  magnitude  and  distribution  of  the  material  properties  throughout 
the  body  is  defined  in  this  work  as  design  synthesis. 

It  is  believed  that  Equations  (9)  and  (10)  have  not  been  previously 
published.  Bert  tl4]  derived  an  equation  similar  to  (10)  wherein  he  reduced 
the  unknown  material  property  coefficients  from  four  to  one.  His  formulation 
is  as  follows: 


[^1  +  2  (c-v)  d4Y 

L  4  r  3  2  2  3  Sr  2  2^.2 

or  dr  r  dr  r  r  Sr  d9 

.  2(c-v)  S3Y  2(c-v+e)  92Y  ,  e  S4Y  ~| 

9  O  '  /  o  '  "T  /  I 


3  *  *.2 
r  Srd9 


S92  '  r4  S94  J 


>1  s3y 


+  dsr  ^i  2^yS?Y_e_aY  2(c-v; 

dr  1/  3  r  _  2  2  Sr  +  2  .  _2 

or  dr  r  r  SrS6 


s! 

r-  99 

i! 

Sr 


_  2  (c-v)-fe  sfi  ]  .  dfs  rai  v  SY  _  v_  S2Yl 

>-3  io2  J  ^2  L-i^2  r  Sr  _2  -,rt2j 


rv  a_  _  e_  il.]  ( 5V  ) 

L  ..2  r  Sr  _2  .  ,2j  \  r/ 


h 

Sr2 

2 

t 

r 


r  992 


(14) 
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with  the  thermal  terms  omitted.  Equation  14)  is  equal  to  Equation  (10) 
with  the  following  identities: 

S  =  a«  n 


e  -  a11/a22 
C  F  a66/2a22 


v  -  a12/a22 


with  S  being  the  dependent  variable  and  e,  c,  and  v  being  fixed  ratios.  It 
is  anticipated  that  most  real  composite  materials  will  exhibit  material 
property  characteristics  as  defined  by  Equations  (15). 

The  solution  of  equations  such  as  (9)  and  (10)  where  there  are  4 
or  more  independent  a^'s,  even  where  these  a^'s  are  assumed  independent  of 
temperature,  is  quite  difficult  but  certainly  not  impossible.  The  simpli¬ 
fying  assumption  made  in  Equations  (15)  leading  to  the  formulation  of  Equation 
(14)  reduces  the  problem  to  only  one  unknown  parameter. 


APPLICATIONS 


Rotationally  Symmetric  Iroblems 


Example  1:  The  pressurized  annular  disk,  internal  pressure:  Consider  a 
pressurized  annular  disk  as  shown  in  Figure  1.  In  the  case  where  the  ring 
is  isotropic  and  homogeneous,  the  stress  distribution  is  as  developed  by 
Lam4  (1852)  and  is  given  by  (c.f.,  Ref.  E 153 ,  p.  60), 


(16) 


where  p  ■=  r/b,  and  k  ■=  b/a.  These  relations  lead  to  the  following  con¬ 
clusions  : 

(1)  |cq  1  >|cTr  |  for  all  p  and  all  k 

(2)  (a0)max  occurs  at  the  inner  boundary  (p  *=  l/k),  and  thus 


^e^max 


(17) 


From  the  stresses  so  generated,  which  for  the  homogeneous  ring,  are  quite 
independent  of  the  material  properties,  it  is  clear  that  the  material  is  not 
being  used  efficiently,  particularly  as  the  thickness  ratio,  k,  increases. 

For  the  homogeneous,  orthotropic  annular  disk  Bienick  et  al.  [ 16 ] 
shows  that  the  stresses  are  given  by 


-(f+1)  .  f-1 

or  *  c^r  +  C2r 


,  -(f+lK  .  f-1 
ctq  *  -c^fr  +  c2fr 


(18) 


1/2 

where  f  *  (a^^/a22^  *  ant*  ci  an<*  c2  are  determined  from  the  boundary  con 

ditions,  in  this  case  ar  *  -P  at  r  *  a  and  a  m  0  at  r  *>  b.  This  work  and 


vork  done  by  Shaffer  [  17 J  show  that  as  k  gets  large  the  limit  of  the  stress 
concentration  approaches  the  orthotropy  ratio,  f,  and  does  so  rapidly.  Both 
Bienick  and  Shaffer  deal  with  the  nonhomogeneous  case,  but  both  assume  a 
special  form  of  nonhomogeneity,  mainly  that 


all  B  all  r* 


a22  *  a22  r 


a12  “  a21  *  a12  r 


where  X  is  ;al  and  a^,  a22»  anc*  ai2  are  constants*  Both  carry  out  optimi¬ 
zation  by  varying  X  and  observing  the  results,  but  no  attempt  was  made  to 
carry  out  design  synthesis  directly. 

Consider  now  that  the  disk  is  composed  of  a  heterogeneous,  ortho¬ 
tropic  material  where  the  material  properties  are  undefined  but  are  functions 
of  the  radius  of  the  disk.  The  equilibrium  equation  for  this  case  reduces  to 


The  compatibility  equation  reduces  to 


—  (ren)  -  e 
dr  v  0'  r 


and  the  stress  strain  relation  for  an  orthotropic  medium  is 


er  "  allar  +  al2a9 


G0  "  a12ar  +  a22t,e 


where  a^,  a^2»  anc*  a22  are  funct*-0Tls  of  r.  Assuming  a  stress  functi 
such  that 


15 


a 


r 


0  dr 


(22) 


and  substituting  Equations  (21)  and  (22)  into  Equation  (20)  and  carrying  out 
the  required  differentiation  results  in 


22 


+  <a'22+^  a^)r  +  i  a'io  A  an)^  -  0 


22 
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2  “11' 


(23) 


where  the  prime  marks  designate  absolute  derivatives  with  respect  to  r. 

Equation  (23)  is  a  single  differential  equation  with  three  independent 
material  coefficients  as  the  operational  parameters.  When  dealing  with  specific 
materials  they  may  be  found  to  be  completely  independent  or  show  some  type  of 
defined  relationship.  As  a  first  step  we  will  assume  that  there  does  exist  a 
simple  ratio  relationship  between  them  that  is  expressable  as 


a22 

*  a22  » 

a '  22 

«=  a'22 

flll 

■  kia22« 

a’n 

“  kla*22 

(24) 

CM 

CO 

"  k2a22’ 

3 ' 12 

*  k2a ' 22  * 

Making  the  appropriate  substitutions  in  Equation  (23)  yields 

I*’  +ir'1']  a,22+|>"  +  rr  -72'']  a22  *°  •  <25> 

Thus  if  Y  is  known  822  Is  fully  defined.  Let  us  suppose  that  for  effective 
material  utilization  it  is  desired  that  cTq  be  constant  throughout  the 
disk;  i.e.,  cjg  ■  Aq.  Integrating  the  second  of  Equations  (22)  gives 

Y  «  A  r  +  B  . 
o  o 

Applying  the  boundary  conditions  that  crr(a)  «*  -  P  and  ar(b)  =  0  yields 


(26) 


with  p  <=  r/b  and  k  *  b/a  as  before.  It  is  interesting  to  compare  Equations 
(27)  with  Equations  (16)  and  (17).  From  Equations  (27)  it  can  be  seen  that 
if  k  is  greater  than  2,  a  becomes  less  than  the  pressure  P,  and  a  which 
equals  P  at  p  «=  1/k,  becomes  the  maximum  normal  stress  in  absolute  value. 
Thus  for  such  a  stress  function  and  geometry  there  is  no  effective  stress 
concentration.  This  is,  of  course,  never  true  for  the  homogeneous  disk. 

With  the  stress  function  now  fully  defined,  Equation  (25)  becomes 


[(1+v..^a.22+[. 


(1~V  kib~[ 

r  r2-1 


‘22 


-  0 


which  has  the  solution 


a22 


(28) 


(29) 


where 


p  «  r/b 


0  ■  -  kx/k2 
k  -k 

_  r  ___i__2 — 1 

s  Ld-k2)k2J  » 


and  for  the  modified  orthotropic  condition  as  defined  by  Equations  (21)  and 
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The  inverse  of  Equation  (29)  or  Eg(r)  is  shown  plotted  in  Figure  2 
for  various  k^'s  while  k ^  was  held  constant  at  0.5.  This  figure  shows  the 
strong  dependency  of  the  modulus  distribution  upon  the  orthotropy  ratio  k^. 

The  ratio  k£  has  a  lesser  effect  as  shown  by  Figure  3.  Here  the  isotropic 
case  (k^  *=  1)  is  shown  plotted  for  five  values  of  k2« 

Example  2.  Pressurized  annular  disk,  external  pressure:  Consider  the 
pressurized  annular  disk  as  shown  in  Figure  1  but  with  the  pressure  acting  on 
the  OD  rather  than  the  ID  as  shown.  If  the  stress  criterion  is  retained;  i.e. 
ctq  ■  constant,  then  the  stress  function  becomes 

t  (i  -  pki 


where  q  ■  external  pressure 
k  ■  b/a 
p  c  r/b  . 

The  compatibility  equation  (25)  becomes 

^2a_l  r(l-ki>  kia”| 

<1+V  '  — J  a'22  +  r”  +  7*]  a22  ■=  0  *  (32) 


Equation  (32)  is  the  same  as  Equation  (28)  except  that  a  replaces  b.  The 
solution  of  Equation  (32)  is 


(33) 


where  0,  §,  and  p  are  defined  in  Equation  (29).  Comparing  Equations  (33) 
and  (29)  we  see  that  they  are  not  the  same,  differing  by  the  quantity  1/k 
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FIGURE  3.  MODULUS  VARIATION  FOR  ISOTROPIC  DISK  WITH  PRESSURIZED 
I.D.,  WITH  <jg  ■  Constant 
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within  the  second  factor.  Thus  the  distribution  of  the  material  parameters 
through  the  thickness  of  the  disk  must  be  different  as  compared  to  the 
internally  pressurized  disk.  This  is  shown  in  Figures  4  and  5.  Also  notice 
that  for  the  internally  pressurized  disk  the  hoop  stress  tends  to  zero  as 
k  tends  to  infinity  while  the  limit  on  for  the  externally  pressurized  disk 
is  q,  the  pressure.  Thus  for  the  externally  pressurized  case  the  material  is 
not  being  as  effectively  utilized  and  perhaps  some  other  criterion  might  more 
suitably  apply. 

If  the  annulus  has  both  internal  and  external  pressure  and  the  same 
stress  criterion  is  applied,  i.e.,  <j  «  constant,  then 


1  *lMp  '  l]  *&[*  *  •,k] 

|j  -  k] 


f 


and  the  material  property  variation  is  given  by 


22 


(34) 


(35) 


Example  3:  The  rotating  disk:  Consider  the  rotating,  uniform  thickness 
annular  disk  as  shown  in  Figure  6.  The  equilibrium  equation  governing  this 
case  is 


tr  <V  * 


,  v  2  2 
+  ■*■  u)  r 
g 


(36) 


The  compatibility  equation  in  terms  of  strain  and  the  stress-strain  relations 
for  an  orthotropic  material  are  given  by  Equations  (20)  and  (21),  respectively. 
Substituting  Equation  (21)  into  (20)  results  in 


0  *  “Vr  +  a12a'r  +  a'22,7e  +  “22" 'e  +  7  1  <a12-aU)<Jr  +  (a22-a12)tJ9 1  • 


(37) 


■  tvJ 
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FIGURE  4.  MODULUS  VARIATION  FOR  ORTHOTROPIC  DISK  WITH  PRESSURIZED 
O.D.,  WITH  a„  ■  Constant 
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FIGURE  5.  MODULUS  VARIATION  FOR  ISOTROPIC  DISK  WITH  PRESSURIZED 
0.  D..  WITH  o„  -  CONSTANT 
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Assume  the  stress  function 


a 


r 


CTe 


r 


2  2 


uu  r 


(38) 


It  is  possible  that  the  material  density,  y,  could  be  a  function  of  the 
radius.  At  this  stage  such  a  condition  introduces  an  unnecessary  complica¬ 
tion  and  will  not  be  considered.  This  being  the  case  the  governing 
compatibility  equation  becomes 


y"(a22)  +  Y-(a'22  +  i  a22)  +  Y  -  -|i) 


Y  2  2  /  ,  ,  ,  a22  a12\  n 

+  i “ r  \a  22  +  3” — r) c 0  • 


(39) 


Again,  for  simplicity,  the  modified  orthotropy  relations  as  defined  by 
Equations  (24)  will  be  adapted  making  Equation  (39)  take  the  form 


0 


2 

0)  r 


(40) 


The  stress  criterion  will  be  the  same  as  before,  i.e.,  CTg  ■= 
constant.  Using  this  criterion  and  operating  on  the  second  of  Equations  (38) 
together  with  the  assumed  boundary  conditions  that  or(a)  e  CTr(b)  *=  0  yields 


(41) 


where 


p  *=  r/b 
k  *  b/a 


v  e  mb  p  tip  velocity 

id  ►  rotational  velocity,  radians 

g  *=  acceleration  of  gravity. 


Here  we  note  that  in  the  limits,  when  k-t  1 


ae 


(42) 


which  is  the  stress 
very  small) 


in  a  rotating  thin  ring,  and  when  k  -*  00  (i.e.,  when  a 


a9 


(43) 


which  is  smaller  than  exists  for  the 
ratio 


max 


isotropic,  homogeneous  case  by  the 


_ 4 _ 

3(3+v) 


(44) 


where  ctd  *  hoop  stress  for  heterogeneous  case 
yl 

ct«  m  hoop  stress  for  isotropic,  homogeneous  case 
o 

v  *  Poisson's  ratio  for  isotropic,  homogeneous  case  . 


Substituting  Equation  (41)  into  Equation  (40)  and  carrying  out  the 
required  differentiation  yields 


»  [~Ar  +  Br  +  D  1 
‘  22  '  Lfr*  +  or2  ♦  HrJ 


22 


(45) 
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D  -  k^(ab)(b+a) 

F  "  k2 

.3  3 

G  -  -(1+k2>  tSt 

H  «  k2(ab)(b+a)  . 

Equation  (45)  is  not  readily  solvable  in  closed  form.  Before  proceeding 
to  the  solution  of  Equation  (43)  by  some  numerical  means,  it  is  convenient 
to  digress  at  this  point  to  discuss  another  stress  criterion  which  may  be 
applied  to  the  implementation  of  design  synthesis  as  defined  in  this 
work.  This  criterion  is  that  through  the  body  of  the  disk,  the  in-plane 
shear  stress,  t,  is  to  be  a  constant.  For  a  body  of  revolution,  acted 
upon  by  symmetric  loads, 


or 


1  ■  <V°r>/2 

x  **  a  -a  «*  Constant  =  C 
0  r  o 


Applying  this  condition  to  Equations  (38)  and  integrating  leads  to 

1  2  2 
-  C  £nr  -  -IW-*-  +  C 
r  o  2g  1 

,  v  2  2 

T*«Clnr  +  C-  4  r  +  C, 

o  o  2  g  1 

4'"  «  —  -  -  Yu2r  , 
r  g  * 


requiring  that  on  the  boundaries; 


o  (b)  “0 
r  o 

a  (a)  =  a. 
r  i 


(46) 


(47) 


leads  to  the  following  relations . 
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2  2 

o  =  =  C  £nr  -  —  +  C, 

r  r  o  2g  1 

2  2  2  2 
v  +  x!i_L_  =  c  (£nr  +  d  ..  m.JL  +  c 


g 


2g 


co  ■  K-i> + 


2  , 
rto  /,  t* 


c.  - - r  £nb  -  o  £na  -  (b  Jlna  -  a  Jlnb) 

i  tn  i  1  »  2s 

a 


These  relations,  as  well  as  those  defined  by  Equation  (41)  (o.  =  constant) 
were  applied  to  Equation  (40)  and  were  solved  numerically.  A  digital 
computer  program,  Determination  Of  Modulus  Variation  JL  (D0M0V1)  ,  was 
structured  for  the  solution  of  these  sets  of  equations.  A  finite  dif¬ 
ference  method  of  solution  as  developed  by  Manson[l8]  was  used  for  the 
calculation  algorithm.  This  program,  which  is  detailed  in  the  Appendices 
was  used  to  solve  several  problems  as  follows. 


•  The  modulus  variation  for  an  orthotropic  disk  with 

a  pressurized  I.D.  for  o^  =  constant  (disk  not  rotating) 

•  The  modulus  variation  for  an  isotropic  disk  with 
pressurized  I.D.  for,  a*  constant  (disk  not  rotating) 

•  The  modulus  variation  for  an  isotropic  disk  with 
pressurized  O.D.  for  ■  constant  (disk  not  rotating) 

•  The  modulus  variation  for  a  rotating  orthotropic  disk 
with  no  edge  loads  for  ol  =  constant 

•  The  modulus  variation  lor  a  rotating  orthotropic  disk 
(Poisson’s  ratio  variable)  with  no  edge  loads  for 

o  =  constant 

U 

•  The  modulus  variation  for  a  rotating  orthotropic  disk 
with  no  edge  load  for  constant  in-plane  shear  stress 


•  The  modulus  variation  for  a  rotating  orthotropic  disk 
(Poisson's  ratio  variable)  with  no  edge  loads  for 
constant  in-plane  shear  stress 


(48) 


•  The  modulus  variation  for  a  rotating  orthotropic 

disc  with  60.0  inch  O.D.,  6.0-inch  I.D.,  turning 

at  2,000  RPM,  stressed  on  the  O.D.  with  an 

uniform  edge  load  of  10,000  psi  for  o  =  constant 

0 

•  The  modulus  variation  for  a  rotating  orthotropic 
disk  with  60.0-inch  O.D.  6.0-inch  I.D.,  turning 
at  A, 000  RPM,  stressed  on  the  O.D.  with  a  uniform 
edge  load  of  12,000  psi  for  Og~  constant. 

The  solutions  for  the  first  four  problems  were  checked  by  use  of  the  closed 
form  solutions.  Equations  (29)  and  (33)  and  are  those  shown  plotted  in 
Figure  2,  3,  A,  and  5.  The  solution  to  the  remaining  problems  are  shown 
plotted  in  Figures  7  through  12.  It  should  be  noted  here  that  when 
boundary  tractions  as  wel_  as  body  loads  are  applied  to  a  rotating  disk, 
the  solution  is  specific  as  regards  the  magnitude  of  the  edge  loads  and 
the  rotational  velocity  of  the  disk.  However,  where  only  one  type  of  load 
is  imposed,  then  the  solution  is  generalized  and  is  independent  of  the 
magnitude  of  the  load.  (Note:  The  stresses  remain  directly  dependent 
on  the  load.)  This  solution  dependence  upon  load  is  shown  very  clearly 
by  comparing  Figures  11  and  12.  Here  the  modulus  variation  is  shown  to 
change  with  the  change  of  the  boundary  tractions  and  the  rotational 
velocity.  All  these  figures  have  been  non-dimensionalized  by  the  expediency 
of  plotting  E/Emax  and  R/B  where 

E  ■  E„,  the  modulus  of  elasticity  in  the  0  direction. 

E  *  the  maximum  E„  calculated  in  the  body  of  the  disk, 

max  0 

R  «  the  radius  of  the  point  in  the  disk  at  which  the 
modulus  is  being  calculated. 

B  =  the  outer  radius  of  the  disk. 

The  computer  output  for  the  curves  plotted  in  Figure  12  is  also  found  here 
in  Appendix  C.  For  this  case,  a.  equals  28,A5A  psi  throughout  the  disk. 

For  an  isotropic,  homogenious,  flat-angular  disk  of  the  same  material 
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FIGURE  7.  MODULUS  VARIATION  FOR  ROTATING  ORTHOTROPIC  DISK  WITH  NO 
EDGE  LOADS,  WITH  a.  -  CONSTANT 
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FIGURE  8.  MODULUS  VARIATION  FOR  ROTATING  ORTHOTROPIC  DISK, 
ORTHOTROPIC  RATIO  OF  2.5  WITH  NO  EDGE  LOADS  AND 
-  CONSTANT 
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FIGURE  10.  MODULUS  VARIATION  FOR  ROTATING  ORTHOTROPIC  DISK 
ORTHOTROPIC  RATIO  OF  2.5  WITH  NO  EDGE  LOADS  AND 
WITH  CONSTANT  IN-PLANE  SHEAR  STRESS 
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FIGURE  11.  MODULUS  VARIATION  FOR  ORTHOTROPIC  ROTATING  ..DISK  STRESSED 
ON  O.D.  RPM  -  2000;  DENSITY  «  0.1  LB/IN.  ;  STRESS  ON 
O.D.  -  10,000  PSI:  O.D.  -  60.0  INCHES:  I.D.  =  6.0  INCHES 
o  -  CONSTANT 
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density  would  equal 


o 


o 


e 

e 


3+v  2,,  2  ,  1-v  2,  ,  0 

—  (b  +  3^  a  )  +  2o 


o  ,2  2 

b  -a 


35,810  +  24,242  =  60,052  psi, 


(49) 


Thus,  by  the  use  of  design  synthesis,  the  maximum  stress  in  such  a  disk 
can  be  reduced  by  a  factor  of  2.11.  Further,  Figure  12  shows  that  this 
can  be  achieved  with  a  modulus  variation  through  the  disk  of  approximately 
2.5  to  1  for  a  material  that  exhibits  and  orthotropic  ratio  k.  of  2.5 
(the  similarity  of  the  two  ratios  is  coincidental  and  bears  no  significance) . 


Non-Svnmetr ic  Problems 

Example  4:  Small  hole  in  an  infinite  plate.  Figure  13  represents  a  small 
hole  in  an  infinite  plate  which  is  subjected  to  a  uniform  tensile  stress,  P, 
in  the  x-direction.  For  the  homogeneous,  isotropic  condition,  the  stress 
distribution  around  the  hole  is  well  known  as  given  by  Timoshenko  [15]  as 


The  maximum  stress  occurs  at  r  «  a,  9  «  (tt/2,3tt/2)  ,  and  is 


‘W  *  ^0)  ,  -  3p  . 

max  0  r-0,e<TT/2 


For  the  homogeneous,  anisotropic  t.  edition  work  by  Green  and  Zerna  [l9  ] , 
Hearmon  [20],  Savin  [21],  Leckhniakii  [12],  and  (as  directly  applied  to 
composites)  by  Greszczuk  [22],  shows  that  the  maximum  stress  at  the  hole  is 
usually  greater  than  for  the  isotropic  case  and  can  reach  values  as  high  as  9P. 


m 


Referring  to  Figure  13,  consider  the  portion  of  the  plate  within 
a  concentric  circle  of  radius  b,  large  in  comparison  with  a.  It  can  safely 
be  assumed  that  the  stresses  at  radius  b  are  effectively  the  same  as  in  a 
plate  without  the  hole  and  can  be  given  by 


(a  )  s»  4  P  (1  +  cos  2  9) 
r  r*b  c 

i  (51) 

(a„a)  «  •  r  P  sin29  . 

rw  r«=b  L 


From  Equation  (50)  it  can  be  seen  that 

(sft)  «  “  P  (1  -  cos  29). 
O  r*b  L 


It  seems  reasonable,  then,  to  choose  a  stress  criterion  for  the  plate 

Og  ■  function  of9=^P(l-  cos  2  9),  (52) 


From  the  second  of  Equations  (11),  with  V  equal  to  zero  the  stress  function 
becomes 

Y  -  J  Pr2  (1  -  cos  2  9)  +  fx(0)r  +  f 2 (9 )  .  (53) 


where  f ^ (9)  and  f2(0)  are  functions  of  9  only.  Applying  the  first  and  third 
of  Equations  (11)  to  Y  results  in 


1  dY  .  I  d  Y 

r  dr  r  2  2 

r  d9 


-  a 


r9 


»  |  P  s  in  2  9 


df2(9) 

d9 


(1  +  cos  2  9) 


l  d2f1(9) 

r  ard92 


d2f2(0) 


d9 


(54) 

fx(0) 


Applying  the  boundary  conditions 


(a) 


r»a 


fare) 


r«a 


0 


yields 
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Choosing  as 


yields 


and  results 


f2<e) 


d  f,(0) 

- V"  +  f,(6) 

d0  1 


Pa 


cos  2  6 


(55) 


-  —  [1  +  3  cos  2  0  ] 

L 


a  particular  solution  to  the  second  of  Equations  (55) 


£^(0)  =  +  C2  cos  ^  3 


(56) 


^(0)  “  “  ^  "  cos  2  ®  ^ 

Y  "  4  {^2  “  2ar)  “  (r  “  a  )2  cos  2  ej- 

CTf  •  j  (l  ~  cos  2  0  +  ^1  -  ~jcos  2  ©J 

(57) 

CT0  *  2  '  “  cos  2  0> 
art  “  -  I  . 

r 
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The  stress  function,  ¥,  as  defined  by  the  first  of  Equations  (57)  was 
applied  to  Equation  (10)  with  the  body  functions  and  temperature  difference 
taken  as  zero.  In  order  to  implement  a  solution,  the  following  relation¬ 
ships  among  the  material  coefficients  t  ’re  assumed. 

all  "  kla22 
a22  “  a22 
a12  "  k2a22 
a66  “  k3a22  ‘ 

A  finite-difference  algorithm  was  structured  to  carry  out  the  solution. 

The  method  of  solution  applied  was  that  usually  referred  to  as  the 
"relaxation  method"  which  is  discussed  in  detail  by  Shaw[23] ,  Hildebrand [25] , 
Allen[26],  and  Richtmyer [27] . 

Due  to  symmetry,  only  one-fourth  of  the  plate  was  modeled.  A 
square  mesh  was  used  and  the  quarter  plate  separated  into  a  square  array 
of  61  x  61  nodes.  It  was  assumed  that  the  circular  hole  has  a  radius  of 
1.5  inches  and  the  mesh  distance,  h,  the  distance  between  nodes,  is  0.25 
inches.  The  width  of  the  quarter  plate  model,  thus  became  15.0  inches. 

This  gave  a  b  to  a  ratio  (see  Figure  13)  of  10:1,  thus  minimizing  the 
far-field  effects  of  the  outer  boundaries  upon  the  stress  around  the  hole. 
This  plate  model  is  shown  in  Figure  14.  Rectangular  coordinates  were 
employed  and  the  differential  equation  that  was  differenced  was  that  shown 
in  Equation  (9).  Equations  (57)  were  converted  to  rectangular  coordinates 
for  input  to  this  program.  In  the  solution  procedure,  it  is  necessary 
to  assume  the  values  of  the  modulus  parameter  at  boundary  nodes. 

The  solutions  were  then  to  proceed  in  an  iterative  manner  until  the 
values  for  a ^  were  determined  at  each  interior  node  in  the  model.  No 
success  was  attained  by  this  method.  The  program  never  was  able  to 
converge  to  a  solution.  In  fact,  a  strongly  divergent  tendency  was  noted 
(i.e.,  each  iteration  on  the  *22'®  at  eac^  no<^e  point  was  markedly  greater 


Outer  Boundary  Points 


Line  of  Symmetry 


Inner  Boundary  Points 


FIGURE  14.  TWO-DIMENSIONAL  FINITE -DIFFERENCE  MODEL 
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on  each  iteration  in  the  relaxation  process  than  the  iteration  which  it 
succeeded,  regardless  of  the  boundary  values  assumed.)  Extensive  investi¬ 
gations  indicated  that  there  wa3  no  error  in  the  solution  process.  It 
appears  that  there  must  be  some  other  governing  relations  which  hav  as 
yet  not  been  expressed. 

As  a  result,  no  solutions  to  the  two-dimensional  problem  in 
design  synthesis  have  as  yet  been  accomplished. 


42 


Limits  and  Other  Considerations 


In  the  preceeding  discussion,  it  was  pointed  out  that  a  solution- 
for  the  two-dimensional  problem  could  not  be  achieved.  This  raised  ques¬ 
tions  as  to  whether  all  governing  equations  have  been  generated.  In  the 
application  of  the  one-dimensional  program,  D0M0V1 ,  it  was  found  that 
under  certain  conditions,  the  only  mathematical  solution  developed  for 
the  modulus  variation  in  a  disk  would  require  that  at  some  points  in  the 
disk  body,  the  value  of  the  modulus  must  be  negative.  Though  this  is 
perfectly  acceptable  from  a  mathematical  viewpoint,  it  obviously  cannot 
be  implemented  from  the  viewpoint  of  real  materials.  To  illustrate  this, 
consider  the  case  of  the  rotating  annular  disk  with  boundary  tractions 
when  subject  to  the  stress  criterion  that  the  in-plane  shear  stresses 
throughout  the  disk  must  be  constant.  From  Equations  (48) 


2  2 


Co(tnr  +  1)  -  *»-£-  +  Cj 


2  2  2 

o  “  C  £nr  -  ■  +  C 


2g 

2g  +tl 


C  -  1  :  |(a  -a  )  +  *”(b2-a2)} 

o  ln  b  l  o  i  2g 

a 

2 

C  b  --■  io.£nb-o  £na  -  -  (b^£na-a2£nb)[ . 

1  £n  -  l  1  0  2g  ' 


(48) 


a 


Take  the  condition  where  the  disk  is  not  rotating;  i.e.,  u>  -  0,  then 


Co  " 


Cl- 


°o“°i 

£n  — 
a 

o.)lnb-a  £na 
1 _ o 

mi 

a 


(59) 


oQ  -  CQ(£nr+l)  + 


o  ■  C  £nr  +  C.  . 
to  1 


(60) 
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Substituting  Equation  (59)  in  (60)  yields 


— — r-  | (o  -o.)(£nr+l)  +  a.£nb-a  £na  i 
in  —  °  1  1  °  > 


Qj.  “  — ^  |(aQ-ai)Jlnr  +  Ojinb-c^inaj  . 
a 

Now,  further  assume  that  oq  -  0,  a  =  -P  .  Then,  Equation  (61)  become 


(61) 


in 


(62) 


in 

a 

Note  that  for  all  values  of  r  <  b,  o  <  -P  in  the  algebraic  sense. 

b  ^  i 

However,  if  —  >  e,  (2.71828),  then  at  r  •  a,  in  —  >  1  and  o  is  negative, 
u  a  0 

Defining  the  displacement  at  r  *  a  as 


U 


U 


U 


(r-a) 

(r-a) 

(r-a) 


(e(6) 

1 

E(6) 


1 

ITeT 


l06'V0r3)  I r-a 

!~T  [l  -  in  £  +  v  in  £] 

( in  — 
a 

I  T  11  -  *n  7  •  (1"v)]| 

(in  —  ’ 

a 


(63) 


From  Equation  (63),  if  in  ■-  >  then  the  term  in  the  brackets  is  negative 
and  U  or  E(0)  must  be  negative.  If  either  is  negative,  then  the  system 
must  do  negative  work,  which  is  not  possible  for  real  materials.  This, 
of  course,  can  be  overcome  in  a  numerical  sense  by  simply  requiring  that 
only  those  solutions  are  acceptable  which  yield  a  material  coefficient 
matrix  that  is  positive  definite  at  each  point  in  the  body;  i.e.,  the 
following  conditions  must  all  be  met: 
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Such  chat  Equation  (65)  becomes 


U  -  II 


Following  the  method  of  Weinstock[28] ,  the  extremization  of  (68)  is 
affected  by  forming  the  integral  1(6)  by  replacing  ¥  in  the  integral 
of  (68)  by 

i  -  'J'(x.y)  +  5n(x,y)  , 


where  ^(x.y)  is  the  actual  extremizing  function  and  n(x,y)  is  an  ar¬ 
bitrary  function  that  is  twice  continuously  differentiable.  Then,  the 
integral  1(6)  is  an  extremum  for  6  -  0,  so  that 

I’(0)  -  0  . 


Writing  f(T  ,  T  ,  V  )  ■  the  integrand  of  Equation  (68),  and  here  the 
xx  yy  xy 

subscripts  refer  to  differentiation  with  respect  to  x  and  y.  Then, 
according  to  (68)  and  using  (69)  to  compute 


>y  3T  3<F 

xx  yy  xy 

36  “  nxx’  36  "  '’yy *  36  “  nxy  * 
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Since  (73)  must  hold  for  arbitrary  n,  it  must  in  particulfr  hold  for  those 

n  for  which  n  =  n  *  n  =  0  on  the  boundary  C.  For  such  n  Equation  (73) 
x  y 

reduces  to  the  well-known  compatibility  equation  as  previously  derived 

(see  Equation  (3))  and  the  resulting  (Equation  (9)).  For  arbitrary  n,  n  ,  n  , 

x  y 

other  than  zero,  the  boundary  relations  must  be  derived  from  the  remaining 
line  integral  portion  of  Equation  (73).  The  first  part  of  the  line  inte¬ 
gral  becomes,  along  C  on  y  =  constant. 


8all  82'}'  a3y 

■  3y  2y2  Sll3y3 


^12^  1  334>  1  3a66  32f 

3y  ax2  (a12  2  a66>  3x23y  2  3x  3x3y 


n[  jV  aVl  iT  aVl 

yLU3y2  123x2J  2La663x3yJ 


The  last  two  terms  in  Equation  (76)  are 


-\[£x]-[\ky]  • 


If  these  strains  are  not  prescribed  as  zero,  then  along  this  portion  of 


the  boundary  n 


0,  and  the  boundary  equation  which  must  be  satis¬ 


fied  is 


3S11  .  a.  ^12  af*  .  ,  ,  1  .  3  V 

3y  3y2  3y3  *  3x2  12  2  ‘66)  3x23y 


+  1  !!«  .  0 

2  3x  3x3y 


and  a  similar  set  of  relations  can  be  written  for  the  second  part  of  the 
line  integral.  These  equations  have,  at  best,  a  very  limited  application 
other  than  showing  that  other  constraints  do  exist  on  the  boundary. 

Their  limited  application  is  due  in  great  part  to  the  choice  of  a  rectan¬ 
gular  coordinate  system.  All  attempts  to  date  to  express  these  relations 
in  terms  of  other  coordinate  systems,  in  particular  those  using  generalized 
normal  and  tangential  components  have  been  unsuccessful.  More  work  must 
be  done  along  these  lines. 
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Material  Considerations 


Finally,  just  a  brief  note  on  the  possibility  of  achieving  a 
material  with  tailorable  properties.  It  is  sufficient  here  to  say  that 
work  by  Adams  and  Tsai[29],  Dimmock  and  Abrahams [30] ,  Hewitt  and  Malherbe[3l] 
Halpin[32] ,  Halpin  and  Pagano[33],  Kohn  and  Krumhansl [34]  Tabaddor [35] , 
Fotinich[36] ,  Wang[37],  and  Fokin  and  Shermergor [38] ,  among  others,  have 
clearly  established  that  material  properties  for  composites  of  various 
types  can  be  established  by  a  knowledge  of  the  known  properties  of  the 
constituent  materials,  their  orientation  and  their  relative  density. 
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SUMMARY 


Mathematical  design  synthesis  has  been  shown  to  be  possible  in 
certain  specific  applications.  The  selection  of  a  design  criterion  in 
the  cases  discussed,  two  dealing  with  stress  distributions,  and  the 
development  of  the  material  property  distribution  within  a  plane  body 
such  that  compatibility  is  satisfied,  appears  to  be  a  rational  basis 
of  design  for  composite  materials.  Difficulty  has  been  encountered  in 
the  solution  of  two-dimensional  problems  employing  this  concept  due  to 
the,  as  yet,  undefined  boundary  restraint  requirements  which  affects 
the  selection  of  the  stress  criterion  to  be  met. 
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APPENDIX  A- | 


INPUT  DATA  FOR  PROGRAM  DOMOV1 

Five  cards  constitute  the  input  data  for  DOMOV1  for  each 
variation  to  be  studied.  These  five  cards  represent  three  (3)  data 
sets.  As  many  runs  may  be  stacked  behind  the  initial  set  as  required. 
The  computer  run  will  terminate  when  an  End  of  File  card  is  read. 


READ  (5,  10)  (ITILU,  1),  I  -  1,  24) 
10  Format  (8A10) 


Any  information  can  be  placed  on  these  three  cards.  It  is 
necessary  to  have  three  cards,  although  any  can  be  left  blank. 

Data  Set  2 

READ  (5,  20)  RO,  RI,  RPM,  DENS,  N,  I BOND,  1ST 
20  Format  (4F  15.0,  3115) 

RO  ■  Outer  Radius,  Inches 

FI  **  Inner  Radius,  Inches 

RPM  ■  Revolution  per  minute  of  Disk 

DENS  »  Material  Density  of  Disk,  lb/in. ^ 

N  ■  Finite  Difference  Increment  and  Printout  Number 
IBOND  ■  Type  of  Boundary  Condition  on  Inner  Radius 

0  -  A  stress  condition 
1  *  A  displacement  condition 

1ST  ■  Type  of  Stress  Criterion  Chosen 


1  ■  Constant  theta  stress 

2  -  Constant  in-Plane  shear  stress 


Data  Set  3 


If  I BOND  Equals  0, 

READ  (5,  40)  SI,  SO,  ORTHO,  PO,  ETHETA 
40  Format  (2F  10.0,  3F  15.0) 

SI  ■  Radial  Stress  on  Inner  Radius 

SO  -  Radial  Stress  on  Outer  Radius 

ORTHO  -  Orthotropic  Ratio,  E^/E^ 

PO  ■  Poisson's  Ratio 

ETHETA  =  Modulus  of  Elasticity  in  Theta 
Direction  at  the  Inner  Radius 

If  I BOND  Equals  1 

READ  (5,  30)  UI,  SO,  ORTHO,  PO,  ETHETA 
30  Format  (2F  10.0,  3F  15.0) 

UI  «  Selected  Radial  Displacement  of  Inner  Radius 


ooooooooooobononooonooopopoo  i  '  o  o  o  o 


Reproduced  (rom 


APPENDIX  B 


LISTING  OF  PROGRAM  D0M0V1 


PROGRAM  00M0V1  ( INRUT, OUTPUT  ,  T  AP  F5  =  I NPUT , TAPE6=0UTPUT >  00M  1 

- C  .  - -  _  OOM  2 

OIMENSION  R(1  M»  ,rTMrT  <10**  )  ,EOC(i:  0)  ,A22(1D2>  ,  CC  (ICC  )  ,00  (1C21  ,  OOM  3 

-  -S . . PHICi;:),PHI°(i:3>  ,PHIPP(lCC),SIGR(10C),SIGT(lCt)  ,U(i:ii)  DOM  A 

t  ,R0°(lw>.)*ITTL(24,i)  00M  5 

_ DIMENSION  39(liii)  .AAll-DKEEliaOJ _ 00“_.  6 

00H  7 

-  -  - - -  .  _  COM  S 

COMMON  R,PHI,PHIP,®HIPP,SIGR,SIGT,U  00“  q 

- - OOM  1C 

OOM  11 

_ 5.J7EA0  .<5ti:».  CITILtI.l).,l=±,2Vl _  OOM  .12 

10  FORMAT  (9Ali)  OOM  13 

- IF  (E  0F  ,5 )  12.,  15  ..... _ _  _  .  OOM  1 A 

15  RFAT  (5,2.)  90, 9 I , R®M,  OSNS ,  N , I ODMO ,  T ST  OOM  15 

.20  FORMAT  (AF15.:,3T5)  .  _ _ _  ..  00“  16 

IF  ( IOONO -1 )  35,25,2*'  OOM  i/ 

_25  READ  (5,3.)  UI.SO.PRTH.O,  PO.EJMET  A _ OOM  1* 

3C  FORMAT  (2F1C.C,3F15. .)  •  OOM  IP 

90  =OUT-P  RAOIUS ,  INCHES  _  _ _  _ _ 00“  21 

91  =imi“R  otQIUS,  INCHES  00“  21 

- 9PM  .  =REVnLUTT ONS  d-o  MINUTE  OF  OISK  _  _  OOM  22 

OENS  =MATE°I AL  DENSITY  OF  OISK,  L0./CN9IC  INCH  DOM  2? 

- 1 - N - .=  F  I  MI  TE  .  C  T  p F E  RE  NC  E  INCREMENT,  NUM9EF. _ _ D0M_-2a 

IDOND  =TYPE  Of  OOUNOARY  CONDITION  ON  INNER  RADIUS  OOM  25 

0  =  A  STRESS  CONDITION  _  _  OOM  26 

1  r  A  DISPLACEMENT  CONDITION  OOM  27 

..  1ST  =  TYPE  OF  STRESS  CPITERICN  CHOSEN.  PRESENTLY  00“  29 

1  *  CONSTANT  THETA  STRESS  OOM  2< 3 

_ _ _ 2.  =.  CONSTANT .  IN- PLANE. SHEAR  STRESS _ _  DOM  33 

UI  =RAOTAL  DISPLACEMENT  OF  INNER  RADIUS  OOM  31 

SI  *9A0TAL  STRESS  ACTINO  ON  INNF^  PAOTUS  00“  32 

SO  rRAOIAL  ST°ESS  ACTTNP.  ON  OUTER  RADIUS  OOM  33 

.  -  ORTHO  ....  =0RTH0T°0oIC  RATIO,  E-SU9-R/E-SU9-THFTA  .  OOM  3A 

p0  spOI SSONS  RATIO  IN  THE  THETA-R  DIRECTION  OOM  35 

- ETHET.A _ =_M00ULUS.0E_,:LASLIC7.LY_IN  J_HE_Thfja_  DIRECTIOM.AT _ 00“_.36 

THE  TNN^R  RAOIUS  00“  37 

.—  ROD  =RATIO  OF  THE  RAOIUS  OVEP.  THE  OUTER  RADIUS  00“  3! 

EOC  =RATIO  OF  THE  MODULUS  OVER  THE  MAXIMUM  MOOULUS  OOM  3® 

.  FOUNO  IN  THT  BODY  0F_  T.HE.OJSK _ ... _  OOM  U'. 

PHT  =THE  STRTSS  FUNCTION  00“  (,1 

_ ?HTP _ fTHE_FTRST  pTpivATIVE_OF_THF__STRfSS  FUNCTION  _  OOM  _  42 

PHIPP  =thf  S-ECONO  DERIVATIVE  CF  THE  STRESS  FUNCTION  "  ‘  OOM  4? 

SIGR  sRADIAL  STRESS  _ _  ..  .  00<  4«* 

SIGT  =  CIRP>U“CF  RENT  I  AL  STRESS  OOM  45 

U  =THE  RADIAL  0 IS°LACE“ENT  _ _  .  ..  00“  46 

ETHET  *THE  CALCULATED  0 IRCUMFFRFNTI AL  MOOULUS  AT  FACH  OOM  47 


B-2 


C 

C 


.A22 


CALCUL AT  IMN  POINT  IN  THE  OISK 
_ =1.  C /ETHLT _ 


GO  TO  45 

-35  READ  (5,4.)  SI,SO,ORTHO,PO,ETHETA. 
40  rO°«AT  (  2r16  .G  ,?F15«  J.) 

45  po  =  (-1 . : » »c>n _ 


PI  =  3. 1415924';  4 

OMEGA*  (ROM*PIl/C3-.UI _ 

GRHO=(OENS)  *(0MEGA**2.:)/(346.4) 

-_A22(  1)  =1 .2/ETHET  A _ 

ETHET(l) =ETHFTA 

_ R(l)  =RJ _ 


OEL=  ( RO-  R  T  )  /  (  N-l ) 

_ .00  5  0  I  =  2,N  _ 

J=l-1 

R(  I) =R ( J  ) +OEI _ 

5C  CONTINUE 
_ RA=_RO/°L _ 


CALL  STRESS  (N, RA ,GRHO ,  SI. SO, UI , O^THO, PO,E THETA, I  POND , 1ST , R I , RO) 
00  55  1=1,  N  .  .......... 

aa  (i»=phip  »n  ♦  (Ro*oHn  m  /  (Rd  n  ♦  (grho*  jrc  i)  **2.  on 

.  OB  (I)  =PHI°P(H  +?HI»(I)/R(I)-.  (ORTHO*  PHI  ( 1) )/ (R(I )  **2.E  >  ♦  ( GRHO*  R 
$i*(3.:-°on 

._cc  a  > =03.aiiLAA.(u _ _ _ : _ 


55  CONTINUE 

_  00(1)=?.: 

EE (1 ) =C«  : 

_ D0.6:  1=2, N _ 

J=I-1 

_ oo  ii  >  =  i  *  :*cc(i)*(p(i)-.p(j))/(i,u. 

FEm  =  l.'-CC(J)M=-(I)-R(J))/2.: 
A22(I)=(cE(T)/C0(I))*(A22(J)) 
ETHET(I)=  1, j/422 (I) 


6C  CONTINUE  _ 

EMAX= 

_ 00.55  .I  =  i,.N _ 


EMAX=AMAX1 (ETHET ( I) ,EHAX) 

65  CONTINUE  ..  _ _ _ 

00  7  0  1  =  1,  N 

......  R00(I)=Rm/RO  _ _ 

EOC(I)=ETHET(I)/phax 

- U(I)  =  A22  ( I )  IjlllJJLLS  LGI.. I J.6P.Q2LS.1.G  2.11)  _L 


70  CONTINUE 

P0=-1 . 0*  °0  .  _  _  _ _ .  .  .  „ 

WRITE  (5,75)  (ITIL(I, 11,1=1, 24) 

...  .75  rop.HAT  (ihi,8ai:/ih  ,(?ai;/ih  ,8Ai:///> _ _ 

ir ( iooNo-1)  «:,r,5,i'»5 

- 8t  WRITE  (6,95)RO,RI,R:»N,5<),si,ETHETA,®0,Oc!THO,OENS,K  _ _ 


OOH 

44 

_ DOM 

49 

oom 

5C 

-  OOH 

51 

OOM 

5? 

_ DOM 

52 

DOM 

54 

—  DON 

55 

OOM 

56 

—  OOM 

57 

OOM 

55 

_ DOM 

59 

OOM 

6: 

...  OOM 

61 

OOM 

62 

...  DOM 

62 

00* 

64 

_ OOM 

65 

OOM 

66 

DOM 

67 

OOM 

6* 

(I  00u 

6* 

OOM 

7: 

_ ocu 

..71 

OOM 

7? 

_  OOM 

73 

DOM 

74 

._  ..  OOM 

75 

OOM 

76 

_ _ DOM 

_  77 

OOM 

7» 

OCM 

79 

OOM 

5C 

_.  DOM 

51 

OOM 

52 

_ DOM 

_  5’ 

OOM 

44 

.  DOM 

85 

OOM 

46 

._  OOM 

5? 

oom 

8? 

_ OOM 

56 

DON 

i: 

OOM 

91 

nou 

R? 

00M 

92 

OOM 

94 

OOM 

95 

B-3 


65  FORMAT  (5X»5  3HnETrok*TR4TTO*l  OF  MOOULUS  V  ARIAT  TON  IN  AN  ANNULAR  OIS  OOM 
SK,//,1X, l.HTNPUT  DATA,/,1Y,1CH********** , //, 5X , 1 5HOU TER  RAOTUS  =  ,  nOM 
SF10.4,2X.15HINNER  OAOTUS  =  ,  Fl  :•  .  4 , 2X  ,6H°PK  =  ,  FI  i  .  2, // ,  5X  *  3C  HR  AO  I  A  TOM 
tl  STRESS,  OUTFR  RAniUS  =  ,  Fl  C  .  2 . 2  X  ♦  3  l'MRA  PI  Al  STRESS,  I  NNEP  RADIUS  DOM 
%s  »Pi:  .2<//t5V,25L'£THFTA  AT  INNER  RAOTUS  =  , CPE  1 5 . 7 , 2X , 1 7W°0T SSONS  DOM 
_ %  RATIO  =  ,FlJ.9,2X,2:H0RTH0Tt>0PIC  RATIO  =  »FlC«4,//»  5X,19H MATERIAL  flO-M 


*  DENSITY  =  ,F1J.5,2X,26NNUMDER  OF  RAOI  AL  °OINTS  =  ,13,///)  OOM 

9C  WRITE  (6,95)  -  -  .  -  - -  -  -  -  D°M 

95  FOPMAT  (  IX,  ll«OUTt>UT  0  AT  A  ,/,  IX  ,  11H*»  •********,//,  7X  ,  6HR  A  01  US  ,  9X  ,  6U  DC* 
SETHETA,  9X,VHSIGMA-P, 9X , 7H3TGUA-T  ,  1*X  ,4HR/P0,  1C  X, 4HE/E0, 1  OX, 7HU-SUD  OOM 
*_ot/, 7*, f,u******,qy,r,H  »•****, 9X,7H»******,8X,7H»******,1CX,4H****,  OOM 

_ S12X,  4H****,1-X ,  7H*»  »*•**,//) _ DO41 

WRITE  (6,1  JO)  (R(I)  ,CTHET  (T)  ,  SIGR  (I)  ,SIGT  (I),  POKI)  ,EOC(I )  ,U(  I  > ,  1  =  1  OOM 

S,N)  .  - -  — - - D°M 

10t  FORMAT  (S  .,P1. .4,2/, 0PE13.6, «X,F7.C, AX,F7.  ’,6X,F«.5,6X ,P9.4,6X,  DOH 

-  ..  SE1  2 . 6)  n0M 

GO  TO  116  OOM 

.10  5  -  WRITE  (  6, 11 J)  P.  0,  RI  ,RnN,SO,UI  ,E  THETA,  PO,  ORTHO,  DENS  ,N - DOM 

11C  FOPM AT  <5Y,52HO::TrPvINATnN  OF  MOOULUS  VARIATION  IN  AN  ANNULAR  OIS  OOM 


$K,//,1X,  1*HINPUT  r'ATA,/,lX,10H*****»«***,//,5X,l5HOUTER  RADIUS  =  ,  DOM 
SF1C.  4, 2X,  15HINNF  R  RAOItlS  =  ,F10«4,2Y*6HRPM  =  ,  Fl  <J .  2  , // ,  5  X ,  2  2  HR  AO  I A  OOM 

_ SL  STRESS,  OUTER  RADIUS  =  ,  F 1 C/ .  2 . 2X  ,  3EHRA  0 1  AL  DISPLACEMENT  ,  INNER  R  OOM 

SAD  IDS  =  ,F10  .7,//,6X,?6HTH-:tA  AT  INNER  RADIUS  =  ,  C  PF1  5 . 7 , 2X  ,  1  7M°0 1  nOM 

_ tSSOMS  RATIO  -_,Fi: .3 ,2X,2:HORTHOTpOPIC-RATIO.  =_,F13. 4,//,sx,l°HMAT  DOM 

SERIAL  DENSITY  =  ,F10«5,2X, 26 HNUMOER  OF  RADIAL  °OINTS  =  ,13,///)  OOM 

....  GO  TO  90  DCH 

115  CONTINUE  00,1 

.  ..  GO  TO  5 _ — - D0M 

12C  CONI INUE  D°M 

_ C/'LL-  _ _ D0‘J| 

ENO  D0M 


96 

97 
9? 
99 

1  j  ' 

101 

102 
103 
10  4 
105 
10  6 
10  7 

10  9 
109 
11C 
111 
112 

113 

114 
1  IE 
Ilf 
117 

1 1  9 
119 

12  0 
121 
122 

123 

124 

125 
12  6 


SU OR OUT  I NE  STRESS  <N,RA , GRHO, SI ,S0,UI,0RTH0,P0,ETHETA,I80ND,IST,RI, 
SRO)  _ _ _ _ 

01  MENS  ION  R(i:a),ETHET  (1C:i.cOC(1C3I,A22(1CO),CC(13CI  ,00(1201, 

»  PHI (1C  J  )  ,£>HIo(l“;  )  ,  PHI  pP(lC:),SIGR  (13C  )  ,SIGT  1101  »  ,U(100) 

- ,R0B(1C:),ITIL(24,1)_- 


DI  HENS  ION  BR(12C>  ,  AAI1  I'll,  EE  110  2  > 


c  _ _ _ 

c 

- -  _  COMMON  R»DHI,PHIP,PHIPP,SIGR,SIGT,U 

c 

_c _ 


STR 

STR 

STR 

STR 

STP 

STR 

STR 

STR 

STR 


1 

2 

3 

4 

5 

6 
7 
6 
9 


C 

C 

c 

.r 

c 

-CL 


C 

C 

C 

c 


c 

c 

c 

c 

c 

_c 

c 


c 


THIS  SUBROUTINE  CALCULATES  THE  STPTSS  FUNCTION  AND 

ITS  FTRST  ANO  SEC0“0  DE RIVA T I VE SFOP  THE  FOLLOWING  CONDITIONS* 

IF  I  BOND  =  D  ANO  1ST  =  1  WE  HAVE  A  STRESS  CONDITION  ON  THE 

-  INNER  RADIUS  AND  A  STRESS  CRITERION  OF  CONSTANT  THrT A 

STRESS 

_IJi.  1P0ND..=  .3.  ANO  1ST _ =  _2.HE._HAVE  A..ST°ESS  CONDITION  ON  THL_ 

INNER  RAOUIS  AND  A  STRESS  CRITERION  OF  CONSTANT  IN-PLANC 

_  ....  SHEAR  STRESS  _ _  _ _ 

IF  190*10  =  1  ANO  1ST  :  1  WE  HAVE  A  OISPLAOEW£  NT  CONDITION  ON 

- THE  INNER  P.ACIUS  WITH.A  STPESS  CONDTION  OF  CONSTANT 

THETA  STR'SS. 

_LH£S£_  ARE  ..THEL_CtlLT_C.0R0IlI.0_NS_£30CPAMMEQ _ 


IF<IR9N0-1)5,35,3’  ..  _ 

5  GO  TO  (1 j , 22 ) *  1ST 
1C  CONTINUE  _ 


C 

C 


- THIS  SECTION  Oc__THE  .SUBROUTINE.  CONT  A I  NS.  THE.  C  ALCUL  A  TI ONS _ 

FOR  A  CONSTANT  CIRCUMFERENTIAL  STRFSS  CONDTION 

AT. a  M  ,  0  /  (RO-Rin*  I  (SO*  RO-SI*  on  ♦  ( (GRHO/3  .u  )  *  CR0**3 ,  r-RI**T .  J  ) ) ) 

_  BP  =  (  (RO*RI)/(RO-R.I))*USIrSO)-.(SRHO/2.0J*(RO**2*C-RI**2.C)  )  _ 

DO  15  I*l, N 

_ PH  I  ( I  >  =  A  J  * R  ( I )  -_( G  >  HO i / •» .  r  )  *_(  ijj  i )  *  *  3 . ;  >  »  bu _ 

PHI^ (I)* A. - (GRHO) • (0(i j **? , j j 

. —  PHIPP(I> * (-2.CJ* (GRHO) *R(I) _ 

SIGR(I)*RMI (II/RCTI 

—  SIGT(I)  =  °HIo(I)4-GRH0*(R(I)**2«0) _ 

15  CONTINUE 

_ GO.  TO  4  j _ _ _ 

2L  CONTINUE 


THIS  SECTION  OF  THF  SUBROUTINE  CONTAINS* THE  CALCULATIONS 
.FOR  A  CONSTANT  IN-PLANE.  SHEAR. STRESS. CRITERION  _ 


- STR 

1C 

STR 

11 

_ STR _ 12 

STR 

12 

STR 

14 

STR 

15 

.  STR 

16 

STR 

17 

_ STR  . 

19 

STR 

19 

_  STR 

2C 

STR 

21 

STR 

22 

STR 

23 

_ STR  _ 

24 

STR 

25 

_ STR 

26 

STR 

27 

.  ..STR 

29 

STR 

29 

_ STP. 

7  ^ 

STR 

31 

STR 

32 

STR 

1 3 

_  STR 

34 

STR 

35 

STR 

STR 

37 

__  STR 

T9 

STR 

39 

_ STR 

41 

STR 

41 

_ STR _ 

42 

STR 

4? 

STR 

44 

STR 

45 

...STR 

46 

STR 

47 

V 


tMl 


Reproduced  from 
best  available  copy. 


033333  J393033  o  3  3333333^333  3|3 

«aJ •  J  il<  J  al  J  J  J  uJ  J >  il  jJ' il  iiol  ~iJ  J  al  'ii  ‘ J  i  J  il  iJ  iJ  'ij  uJUJ 1 

<\j  -«  t  m  n  d  -«  vj  -o  -r  ^  ^  i  o  -d  -O  o  -•  4-  >  \i  o  "O  ♦  r» 

K-n,-*  >  'me  -*  m  ^  n  o  .n  o  o  --  *  o  \j  o  o  -« n  oko 

M  />  -*  O  3  >  ♦  M  "V*-  »>-*-•  -*^V  O-  Z>  *  >  /I  -•*-  *  -•  X)  /»  nio 

>J  OlMtOlOintVI^  ,».>.>. >.>«*.*!>.  -3- -O  .3  -•  -M.M  Oif^- 

♦  iro  m-VJ  M  M  M«  VJ  M.-*  -*-«-«  -*  \l  ■'J.'VJ.M'M.M  M  M  M 


O  O  300010000  o 


OO  O  O  OOlOOOOiOOlO 


o  -#■  oo  m  co  a?  o  'M1  to  qd  #  or»  o  nj  m  -o  f** 

n  o  />  #  ni\J  -•  >'*-  in ^  t>  ^  o  o  \j 

n  m  -oir w  x,  o  o  h  ,m  n  ♦  ^  m  «x>  r*-:o  co  o  o  o  ^  rvj  rvj  r>  .4  •*  in 
-•  efvirvjrg  vjrvjrvjoj  rvj rvj  r\j oj  mm  mr>mr»  nf»inn 


zoo 
*  o 

M  -< 

•OH 

SI  I  —  7»l 

UJ  <  »- 

x  x  z 

—  -« 

y)  /)i  3 

!  z  a. 

J  O 

<  n  j 

•■*»/>< 
.3  M 

<03 
X  X.  < 

X 


•  i  ••  •  •  •!  •  •!  •  •  •!  •  •  •■  •  •)  •  •*  •  •»  *i  •!  *•  •  I  ••  •{  *i  •  •• 

♦  ■*****<*■•* 

in  nin^n.o.AJi  j)inin  j>ninj)j)j>  on  nix>  .n  .n  n  ji  in  n 

4  t44444tt|4  4444  4.44  4  4  t  #|4  4  4  4  4  4  4 

0.00  o  o  oooon>  o  o  ooolooo  ooo  oo  o  oo  oo 

M|M  JV  M*M  IM  ■MlM  V  \l  MM  \l  'MMIM.M.M  V  \l  rvj\l  "M  NJ  M  M  fM  fM 


Or  •  •! 

:♦  r>l 

k> 


C  <H 

5  <i; 

H 


3  H 

< 

X  in 

D 

X  -« 

U  3 

—  < 

d  r 

o 

X 

•  J 

o  ^ 

A  Z 

il 
X 


J  4 

*  — 

3  X 
<  *- 
x  a 


4  J  3  ;  3  Jj3  i  ;  3  3  3  3  3  ;  il  ;  WSJ 


*  • 

3  * 
I  -«  • 

I  3  * 

I  *  • 


it  o  no* 


I 


3300000  3  D  D  3  0  0  0  3  3  0  3  0  0  0  0  3  00000-30000  000000000000  30 

•  •  •  •  1  •  •  •  •  I  •  I  •  •  •  I  I  •  I  •  I  II  II  I'll# . lit . 

J  J  J  J  iJ  iJ  J  J  J  '4J  Joi'ii  'U  J  iJ  J  J  Xj  jJ  il  xJ  xJ  xJ  il  J'  xJ  ii  xJ  J  xi  U  xJ  J  xi  U  xJ  J  JxJ  xJxJxJJxJxJ  J 

3>  •*>  \»  o  *  •'-*'>  '■n  o  m  **vx>  -*o-<j>  ^  *•  o  rwi  ?  ^ 

>  o  oi#  m  ;>  «o  *m  >*■  -*  #/)  x>  n  ’'no  o  >  m^.'vjo  o  — •  #  r>  /v*o  >•*■>-«-  o  oio  m  o  >  o  -n  **-  3  -<  o  o 
ni-<Bo»-tT  on  3 '  »■-*  —  n  o  o  \J  o  -n  >  t  7>  t  >  iomo  3»  i-  »  nano  o  -<  n  />  o 
•  3  m  -n<  »•■»••  n  01— ••»-  oi>03-«  m.m  m  »•■*.«  ft  oioi'o'^oioixx  oo  —  —  —  m  m  u  innin»i«  »i* 
n f  1*  1  * •  «#■ *  •*■  «>#•  t •  <#•# *# »*■  *  *  #  •  />  n  jv/>  n  n  ^  n  n«f>  n*o  n  rt  n /> 


0  3  3  0  0 

I  I  •  9  • 

xJ  |J  xi  xj  xJ 

>  n  o  -«  o 
t  t»  ,m 
o  >  -•  Mi'n 

r  <*  n  nn 
n  si  n  ft  n 


*'  *I*T«  *  ■#  •# 

t  ♦  ♦  ♦ 

OOOOOO 
M  M  M  M  >1 


•  •  • 1  •<  *1  •'  • 

♦  iti-f 

II  SI  SI  Si  Si  SI  S!  SI  SI 

*  *  ♦  *♦*)#.*# 

ooiooooioot) 

V  M  V  >iM  M  M 


♦  ,r|* 

/I  SI  SI  SIVSI  SI  n  Si  SI  SI 

♦  t  •#  *1*  ♦  f  ^  ^  ^ 

000  0)0  t>  t  o  o  o 
vm^vimmv  mm  m 


•  #'*’*«*'  <*)*■  t>*  *  t»f  •  If* 

jw>ji  n  x>  ji  .n  /1  xt  .n  x>  x> 
*  ^  -r  -r  *  »  »  >r  •*(* 
0000  olo  DO  O  O  D to 
lUl.-M  "4  M  M  M  UiMiM  MIMIM 


UxiU*i  jJiJJU’JUkii'uJ'l  i  J 


j  n  \»  >  « 
>  x  n  o 


1 

1 


U« 155<*^<»E^u| 


MilSsOIN*  K A 1  IC  tUUALS  li 


OOOl  00000003^3; 

®  •  •  ■  I  •  i  I  I  I  I  I  I  It 

Ul  Ul'JIJ  oj  i  Ju|,u).jJ:  J  Uii  J  xli 
Mi)'VI-f>NN.O  d>  *  O  XI  ' 

N-rn.x>  .3>-o  >m  a  o  Aio  ***<a  > 

^  n'O-.ODO  »  . 
AJ  -  v+r.*.  A  ,— •  3-t I>ir I  J  ’ 

♦  ♦!/»  Aior—*: 


•  rv  o  o 

m  i 

I  aJ  aJ  jI 

a  >  x> 

A  *-  -<• 

•n  o  A 

■  A.O 
*vi'  aj  aj 


30  O 
III 
_lJ  jJ  x) 
>  O  \J 
oro.-* 
ON-« 
MO!« 
>iO*M- 
aj  nrni 


ooooo 


^  ®  <nl®  ®  -•*.  ,o  air^^iA  x  -•  ONOio-*o^l>r\iNNosU 

^  ^  M  ON  V  A  9  A tA  9iO  ^  ^>  00  ^IC 

9  f\J  9  *AI  lO  N  5  o  0>  >  9  K  J)  #  ■W0>NJ)l*>«4099J)n*4 
><><><*  a  >  O  O  C*  X)  9  999  r-NNNN 


oo  0.0  o  O  oo 


o^  COlf- 
D  3  -*  M 

o.-k.  x»  r-  . 

O-O-^IM. 


o  o  o  o  o 


oo  oo  o 


9  M.O  J>  i/>  ■•OOJI-^OOCO*-  ollfl 

•  m  n  ■»  ji  o  n.  o I >  o  o-«  AJ  "o  t 

l>lO.J>»  Tl\J  h;  jo  >  .  O  •«.  •  O  l/V 

►  \j.  'll* ■  n  oi^. o  >.  j  o-«  xi.niti 

I  AJl  M  \J  AJ  AJ  AJ  AJ  |  AJ  n  'll  1  -*ip 

►  •)•  ••  •  •  •  I  •  •  •  •  ■  •  •[• 

►  ■»  »|3*  ®-»  .3  33  33I3  .^3 


•  •'  •  •  •  *'  •'•'?  •'  •!  •  •'••<•  •!•■•!  .1  •  ..  .1  ..  •  >1  .  .- 

♦  ,-#■  *3  (^'  3  i-3  .«#  13  3  3  ‘3i3  313  3I3  3  13  .3<3  .  3  I3  3  33,3  3I3 

■Q  ■O  -Q^  -J*  -*1  -0  -C  -^  -Q  -ft  si  si  si  ji|x>  si  si  si  si  jaIja  x>  jixi  si  si-si 

**  J  3  I  3  j  ^  j  j  #1  J  J  -J  J  J,  J  -J  *  J  J  r  J  A  A  A  A  A I  A 

J  y  ®!n  8  °  o  o  o  i  o  o  o  o  o  o  o  o  i  o  o  o  o  jo  o  oo  o  olo 

MjAJ  AJ  V  rv»  (V  AJAJ  AJ  AJ  AJ,  \J  :<J  ,AJ  ,AJ|  AJ  AJ  AJ  AJ  AJ  AJ  AJ  \J  AJ  AJ  !.AJ  JAJiAJ 


19  n  f  f  ^  Di’f  n  ^  J)  3  A  f—  f—  ti  ♦  v«  j|  m  « 

1  ^  ^  ^  O  9  J)  <\l  M  **■*  *\i  0*1  ^  sOA  «V  OX  N  ^  >  R 

iT!  T?1:T  *  °  &'**■  *>  \J.r-  j  i  n  o  ojin  ji  m  *\o 

♦(A'^-O  *  t  JIlA  O  O  O^^INN.  ^  D  0»‘0 


O  O  -j#  O  3  0 
♦  ♦♦♦♦♦< 
ii  iJ  jl  -I  j  J„ 
'00.00  » 

3  *n  o,^  «3»  •  a  r 
•O'-  3.|»  *j 

•*  o  nj|t  >  t>  > 

»A  >  |\J  *  0  ' 

o  o  oK'>n> 


■»|  3  3  0  3,  3 


WN  l\  J>  9 

3  g  A  U  3-> 

t  *•  >  n 
■»  w  n  o  i  «) 

•  •  •  •  •  • 
"»;»»  C)  O  f  t 


^ ^  %-s-K  r-  ^  ^  K  ^  ^  ^  t-  ' 

J'iU  blli^j^MJlJ  J  J  U  U  U  U  J  ij  J 

f  >  I  1  ^  I O  O  I  9  ^  ^  )  \i  >  /  f  ; 

1  ■>  V  J  JI  JN-yJ  JJ  JJ  J  >A>AJ>J 
l»»  II  .’I  —  T  ini'-  IM  «  >•)  ^  u  I-  J|  -  J  J|  j 

r  V  -I  VI  -»^  Ai  JU  VI  —  AJI  U  UJ>  I 

»"O.Af|ll-i>ll'/l  t  vl  —  >0' 

V*  '  - - |»  A-r-  3  0  0  |J  J  30. A  I 


3  3  3  3  3|  3  J  O  3  3  3  3  o  > 


'  tiln'i'  jn  ji  m  r  4  »  -  .j./i  S  .  >j-  a, 

,  “I*  f,  >M/UI|>'AOJi-f|'  J  \J  ll  0 
'UJ1VJ  J^J  -  JVOO  >o.vg  }'  J  -,  ■ 

T,  .'I  .’  'i  °.|  :Tl  *  “  '  j .  ;«■;» 

n  n  n  o  o  o  t.,"  —  - -  o lo  j.  >  >  > 


I 


I  -4—  O  o  ODD  9  3  0  0  303  3  3  000000000033 

30000  OOOO  000003  3  300000000000  3  >003  OO  000003000  30  003003  OO 


Lii  J  ii  J  UJ'  J  li 
nMAyo-o©- 
A  A  >  ^  9  y . 

y ^  n  ^  tio  -« 
i  3**  A  AM.  -•  -* 
o-«  A*A  r* ?><«-• 

•  y  y.y.y  •*■»  yu> 


xl  Xl  Xl  P  < 

♦  "V 

o  ©  O  *M  * 
-«♦  M'n  * 
-*  -«  m  ni 
aia 

A  A  A  A  . 


I'Xl  xl  xJ 
•  ♦  y  s-.  - 
» a  o  y 
^  \i  t*  : 
©!-•  A». 
IP  fliO  ’ 
>  O  P'f**  f 


I  xl  Pul  xl 

i  p-vj  a  -n- 

►  O 

i  ^  y  a  n  * 

>  a 

I  rw.^iiNJ 
.  fx. -—•.©;  0)  : 


I  li  I  xl  xl  xl  xl  U  xl  xl  xl  xl  xl  xl  xl  xl  xl  xl  >xl  xl  xl  ui  Jul 

‘  a  a  mo  o  ♦  o  ♦  o  ■n,r>  t  p  tn  ©  tx.  >  o  j»p,n 

■  p-n.M  o  0'*i  a  m  p  ^</>  o  p  t  p  hp  p 

■  OMO  OPP>‘M^PP  >  o 

I  -H  yt©IO*«^  niOlO'H'nm  X..9  \J  tIOlDfH 

>10  0  3  3— *-•  H  h  V  M-M  M  A  At  AiAifi  y*ytf.y  tA 


xl  Ui  xl  Xl 
3  ©  fiM 

m  y  At© 

^  H  H  3 

At  a 


A.  A  x-ov. 
A  A  AjA 


9  0  9  0 

O  O 

OiA  A  © 

o  0000 


00000  30  000  o 


o«h  >  m  onHmnP  m  *y  a  ouya  ©  y  nontMco  y  rvi  — «  r>i©  p  n  h®  p  p  od  piaj  np  m  n^mp  h®  p  pip 

MP.OON^ilvlotiO  ©'O  >-OH  n.P  OMP  OpHO'M  ©  ♦  3^-#  0>'N'0.p'l  y  .  A  iA  •  A  •  A I A  •  A  •  y  ♦ -lA  OK 

o©  r-  -v  y  ro  rvi  — •  o  x>  j  p  p  ^  rAM-<-«:o>3a)®  r»  fx  p  p  p  p  Pi-t  y-  y(A(AfA(0(A.\JMMr\JMM(\ir\iMMrvM\JM<\JM  micm 
kj©  ppp  p  Plop  pp  ppppppp  «PP  ^  t  #  -t  <t  #  y  t  y  y  #  ■#  y  t  t  ■#  »  t  *  tt  y  -y  y  y  y  y  iy  y  *y  y  ty  ♦If 


O  OO  00  00  0900 


0  0,0  0|00|0  OOOO 


OOOO  OIO  OOOO 


O  3  0  OOOjOO.O  OO  0  3 


tflf  AMHOW'ONfPP'tnWHOypNPPP#  (AM  *-•  O  O  ©  rx!©  pp  <t  A  PN  pipp  ^  (A  fMHo^pNKOPP 

A©  r»  ©  O  oantyrn  IPCIS  310  OOHlMntPONpTk  O  C*-«  M  A  y  A  ©  rx.  ©  iO  O  O  -*  f\J  A  y  .fl  O*^"©  OIOO-H.AJ'A  y  A 

y  a«  v-*  >o  >lo.^'on  .♦•oim.h  >  o  >«ahom  »'*»»v.h  ooo*«^»0'a  y  a»\j  -•■  oo>  ©ix.  ©,/>  ytA.  m  -•  31  3  >•©  ^-oip  ♦> 

AOI'-'OIM  3  3I-*  Vi.  At  ft  A  Of  Ot*>  3*  PHig.  •»itin-ON<D'>0  3-<  VIIA»  tiPOl^'  Dt  >♦  3*  3-h  M»  At  *  I  A  -  3rx»  Ol  >1  3- 3--h  VJ  At  y  I A 
(A  A  A  A  “A.  t  f|  *  *1  *  ♦  fflffif  A  A  AIA  A  A  A  A  A|A  A  ©  ©  PPlP  O  P  P  ©  OO  x_  —  fx.-— x_  ^  —  x_  x_ ,  D  jj  0  ©  ©  ©  © 
•!••••  ••••!•  ••  •;  •  •••••  ••••••)••  ••  ••.••••!•  *j  •  •••  ••!••  •  •  «•  •!•  •.••)••• 

03  0.3  O  3  O  30,3  O  O  Oj  3  0  3  0|  O  3  3  OO  OO  9  0  OO  O  O  OlO  0  99.0  Pi  O  0.0  O  O  39  OO  0  09  990  9(0  O  9 


4  •  ••  •'  •  •■••■  4  •  •  *1  •  *1  •!  *1  •  »i  •  •  ■  •  i  ••  •  •  •  ••  •■  •  ••  .1  •  •>  •■  •  •  •  •!  •  •  • 

»)•*  *•*•*<*■'■*  *•*)*•■*  i  *  i+ 1  r -*i*  •  *i* -r  !*•■*}«■•+ !■*•  *ir  *]  +  >*>*■•*<*■  t\-ri*it  r  r  *]*  *.*•**■'*[*• 

stn  s.r>minmi.nmm.n  r><  n  si  n  ^  n  s>\n  si  s>  s\  r>  sti  n  n  n  n  s>  n  s>  n  n  n  n  s> '  n  st  n  n  m  ui  in 

*  #  ♦  *  •#!  *+■*■**+{ +*■»** -r  ***■  +  *  *!»♦♦  t  *  *\-r  ■*  t  **********  t  +  ■*  *  »♦■.»♦**.**(> 

ob  a  a  o.b  x>\x>  x>;b  bo  ui  o  so  oooia  oaa>aaio  a  as  a  ala  as  ssa'ooaas  a  a  o  ass  so  o  os  oslo 

■VI  V  M-M  SMfM|mrvl|  S  M  \J  fVI.  "VI  *VI  M  IV.VJMM  M  M  M  tVI  MI'VI  .M'S  'VI.'V  rvli'M  tM  rvi  Mi. VI  M  M  M  M  M  M  M.  M  M  M  M  \J  Ml  \i  rvi  MMM  M|M 


ins  «.\i»ji.\i  osU»  a 
ms  a  st  Mj'moj- 


am  *»  ji  —•  mi  £•  o  aj>T»  vils  si  u  mv  ii  avh  r  ^  rt  n  vi!s  oms  u  nla  "•  r  v  t>  ah 

m  —  Jio’iij'iiHAsiiK  =•  "s  s  fltua  t  :  x  toms  »■  an  cjii»nN  iw  o» 


’~-.o  ■>  a  —  -a!"in  iin  *  t  i.  n  hvim  m  «  — ■-<  aj.j.  s~-  milt  ,-t  rvi  —  -a  s  t  "tu  u  ssin>S|—  as"S  ‘I  — 

s's  o  s  >  r  >I7<  >•>  a'IS'  >  t>  s  >  >.s  »  >  a>  t>'S  s  ss  a  s|s  s  s  ss~»p'~»"»"»"-'-|'-o-ssosos>nnmrn 


'  3  P  3 
♦  ♦  ♦ 
xi  J  u 
I'  J  3 

II  ©  -H 

^  r  o 
*  r  ju  jj 
M  —  -* 

A  *A  A 


<  9  J  9 
♦  ♦  i  ♦ 
*1*1x1 


ii  «  ^  M)  s  «  g  i  vj 
^  r  3  xi  ;  p  w  p  j)  p 

I*  P  P  P  *  O  '  V!  T  P 

m—  ^  s  o  n  t 

■A  A  A  t  t  tilt  t  T 

•  •••  •••!•  •• 
3  *  3.  >  y  3  3  0  0  3 


3  3  ^03933^9. i  3  333333  333^3 

■xl  tJ  Pi  -J  xl  J  U  u  xl  P  xl  J  U  P*iil  9  'll  xl  <1  OIJ 

f  -*|  .i  g  o  t  i  a  o  -  g  m  m  ')  'i  4  m  ^  1  o  a 

♦  M  f|  O  O  y  A  -•  si  JO  o  ♦  4  A  ©  3  y  1*  .g  it 

j>  p  y  a  ©  m  —  a  ph  ji  ji  n.  .v,  j  Jim  p  'I  -  h  y|p 

i  P  l\l  Jl  l»  A  ^  p  -  u  (P  p  .yl  i  p  A  3  ■*•.  -T  — 

A  y  *a!  a  m  —  —  yp>>o  o  o  o  A  nl.1) 

»rr|TTt»»j'|>in->s  sis  s  s  •>  s  sis 


^  X.  —lx.  -X  X  X  X  ^.X.  fxlx.  XIXX  X  X.  U.  —  —  )x.  s 

3  3  3*  3  3  3  9  3  3  J  .J  3  3  3  9  »  9  3  3  3  9'  O 

♦  ♦♦  ♦♦  ♦♦♦  ♦♦!♦♦♦  ♦♦ 
UPxl9PJUU.il/iilPU  PPPP.  P  .JU  PUI 
-A  J  D  -  -  \J-«0^0|y*‘  O  —  A  0|9  t  f  —  I J  ' 

©4  JM*-D',-A.AO©,MT*  >J»  JMA^I^  r  VJ 

sT  VI  VI  *1  P  J  P  #  A  A  II  O  A  9  U  U>  3  A  U*fl 

f  ^  p  m  -i  X  u  x  ^  j|  r a  **I  a  m  M  vi *  vi  '*i  A  '•i  y  : 

*  *  *  *  *  l  *  *i'n'n^-*i|*v*t,»vn*vnl*vn  AAi*i 


y  y  y  y  #|t  A'aa.aa'aaaaaa'aaaaia 
AAAAA^AAAA  At  A  A  A  A  A  A I A  A  A  A  (A 


3  3,3  3  3  003333  33  33  33  3  0  3  >3  3  3  33330 


K/1,1  3SP.V7P  3  f-  Ji  IVI  >  O  yl—  l|p  AljN  A  MO 

H|  X  JVI  PPf.lAP  >|Hf  ^  J  M  JI  ©  3. A  Ui  >  -*  y  *«-.  j  vi  A  ©  I  J 

ionM  y-<>3  a  ooA'j  j'xy-^^o  a  >  o*n  j  3*x*  ^ 

:  •  ;  'j  «  »i  ?  ;  •  j  ^  ^  ^  j  ^  r;  » 


3  0  3  0  3  3  3  0  3  313 


O  ♦  —  ©  .n  A  *SPM>jO  i  -  P  P  A|  -J  fx  1/,  (\|  ,0  ©  y 


J.vil  viMtAAAAjyy  y  a  a  a  n  o  o  o  o^l-tx..©  ©  ©  ©I 

-  -  H-  -1  H  -•  H  -  -A-  -A-  -A 


a  j  y  x-  3  'Vi  n  ©  .i  A  p  o  -  y  x.  x  M  n  ©' 3  a  3 
3A3©Al^Ax.^  -1  >1  A  A  3  C  A  J1 

A  4>  »  >*  O  U  OA  II  n  v  ••  A|0 


oo>.  3  3!3-#-<-i-*\j'ViMvJAAA|y  yyylA  aIa 
»  -^1  VI  VJi  M  M  M  (M  M  M  M  .VI,  VI  >V«  VI  'vlj  VI  .U  M  M  M  ,VI|  VJ 


J-ii  l Lit.* 


Ha)  1C  £(JUALS  V 


•030' 
♦  ♦  • 

U  oi  iJ  • 
t 

0  O  . 

o  o  - 
»  m  o:  3 

i"-  >IM  ■ 
O  o  p*-  r 


3330>0  0303S; 

*  ♦  ♦  ♦♦  ♦•#  ••  •  4 

'xlxl  J  J  Jij  ij  u  J  j  iJ  A 

_  >  o  ■n  N'-*  ♦  >  -o  o  m  j 

\J  m  n  m  -on 

■no'Oim-i  a  t>io  *!>,>,- 

d  -om'o  o  3\j*  o  r 

ij  Mn'-nnH  iit  ».»  tit 

•  •  •••.••••#•, 

3  00  3  0  3  3  3  0  3  ol  ; 


n.'£  ;?*■  ^  *'■*  —  ^  »  —  *>'■£>  <n  -« o>rw  mU 

2  “?  ——r-'—  __  ^ 

ooooooooooooooooooo 
•••••••••••  ••  •••••!• 

0000  0  0,000  30  O' o  O  30  O  3k3 


^  ®  >°'-n  -n‘*  r>'M  -Js»  O  ■>.  o  /»[« 

* °  ^'2‘n  '»-•  »|  >>'t>  '.a  m». 

3 " •  vi-hi. rin -oi  —-o  >•  >-3  -« \i ,  <>itui 
"*  ---  r-,n  D  D  O  O  D^> 

•  *  •  •  •  •'••••!•  •••... 

3  3  33  33;3  30  310  0  33  3  3  3 


ZZ  ZtZ Z  ZZ Z Z\Z 


n  rt'.i  run  n|n  n  n  r»  n 

>  ->  3  3  3  3  3j  J  3  3  3  O 
l|Ai  J  U  J  J  liLti  J  iJ  J 

*;  *  ^  ^  v  ^  e*"*  -n  -n  n  j)  , 

' !  U  —  —  Ji  Tit)  J  t  x-  —  I 

1}  D  M  ^  ^  U  I  9  ^  p  tj  i 

It  •  u  o  n—  mo  u  t  n  g 

*  ^  *  o -r  s-'j ^  - 

3  >,  »  5  >  >  3  3  j  I  I 


."*■  -Cl'*  fv»t»*lo  »■'—  u  r>r3  ji  \  *  -I, 

n  3  0  n.v,  C C? J;  Z  Cl ?  iJ  %  o  7,  "Z  V  “  / 

.  .; .  .j . .  y  i  °.  :  « ; 7;° 

--.  V3.«.M  Vi  Vi  VIj-Vi.M  .Vi,*.*  -u  Vi  Vi  Vlij^lv  'J 


